Using in silico models to simulate dual perturbation experiments: procedure development and interpretation of outcomes.

BackgroundA growing number of realistic in silico models of metabolic functions are being formulated and can serve as 'dry lab' platforms to prototype and simulate experiments before they are performed. For example, dual perturbation experiments that vary both genetic and environmental parameters can readily be simulated in silico. Genetic and environmental perturbations were applied to a cell-scale model of the human erythrocyte and subsequently investigated.ResultsThe resulting steady state fluxes and concentrations, as well as dynamic responses to the perturbations were analyzed, yielding two important conclusions: 1) that transporters are informative about the internal states (fluxes and concentrations) of a cell and, 2) that genetic variations can disrupt the natural sequence of dynamic interactions between network components. The former arises from adjustments in energy and redox states, while the latter is a result of shifting time scales in aggregate pool formation of metabolites. These two concepts are illustrated for glucose-6 phosphate dehydrogenase (G6PD) and pyruvate kinase (PK) in the human red blood cell.ConclusionDual perturbation experiments in silico are much more informative for the characterization of functional states than single perturbations. Predictions from an experimentally validated cellular model of metabolism indicate that the measurement of cofactor precursor transport rates can inform the internal state of the cell when the external demands are altered or a causal genetic variation is introduced. Finally, genetic mutations that alter the clinical phenotype may also disrupt the 'natural' time scale hierarchy of interactions in the network

What do cells actually want?

Genome-scale models require an objective function representing what an organism strives for. A method has been developed to infer this fundamental biological function from data.Please see related Research article: www.dx.doi.org/10.1186/s13059-016-0968-2

Scalable computation of intracellular metabolite concentrations

Current mathematical frameworks for predicting the flux state and macromolecular composition of the cell do not rely on thermodynamic constraints to determine the spontaneous direction of reactions. These predictions may be biologically infeasible as a result. Imposing thermodynamic constraints requires accurate estimations of intracellular metabolite concentrations. These concentrations are constrained within physiologically possible ranges to enable an organism to grow in extreme conditions and adapt to its environment. Here, we introduce tractable computational techniques to characterize intracellular metabolite concentrations within a constraint-based modeling framework. This model provides a feasible concentration set, which can generally be nonconvex and disconnected. We examine three approaches based on polynomial optimization, random sampling, and global optimization. We leverage the sparsity and algebraic structure of the underlying biophysical models to enhance the computational efficiency of these techniques. We then compare their performance in two case studies, showing that the global-optimization formulation exhibits more desirable scaling properties than the random-sampling and polynomial-optimization formulation, and, thus, is a promising candidate for handling large-scale metabolic networks

Metabolic flux balance analysis and the in silico analysis of Escherichia coli K-12 gene deletions

BACKGROUND: Genome sequencing and bioinformatics are producing detailed lists of the molecular components contained in many prokaryotic organisms. From this 'parts catalogue' of a microbial cell, in silico representations of integrated metabolic functions can be constructed and analyzed using flux balance analysis (FBA). FBA is particularly well-suited to study metabolic networks based on genomic, biochemical, and strain specific information. RESULTS: Herein, we have utilized FBA to interpret and analyze the metabolic capabilities of Escherichia coli. We have computationally mapped the metabolic capabilities of E. coli using FBA and examined the optimal utilization of the E. coli metabolic pathways as a function of environmental variables. We have used an in silico analysis to identify seven gene products of central metabolism (glycolysis, pentose phosphate pathway, TCA cycle, electron transport system) essential for aerobic growth of E. coli on glucose minimal media, and 15 gene products essential for anaerobic growth on glucose minimal media. The in silico tpi(-), zwf, and pta(-) mutant strains were examined in more detail by mapping the capabilities of these in silico isogenic strains. CONCLUSIONS: We found that computational models of E. coli metabolism based on physicochemical constraints can be used to interpret mutant behavior. These in silica results lead to a further understanding of the complex genotype-phenotype relation. Supplementary information

Chlorella vulgaris (Chlorellaceae) does not secrete autoinhibitors at high cell densities

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/141950/1/ajb211559.pd

Systems Biology and Pangenome of Salmonella O-Antigens.

O-antigens are glycopolymers in lipopolysaccharides expressed on the cell surface of Gram-negative bacteria. Variability in the O-antigen structure constitutes the basis for the establishment of the serotyping schema. We pursued a two-pronged approach to define the basis for O-antigen structural diversity. First, we developed a bottom-up systems biology approach to O-antigen metabolism by building a reconstruction of Salmonella O-antigen biosynthesis and used it to (i) update 410 existing Salmonella strain-specific metabolic models, (ii) predict a strain's serogroup and its O-antigen glycan synthesis capability (yielding 98% agreement with experimental data), and (iii) extend our workflow to more than 1,400 Gram-negative strains. Second, we used a top-down pangenome analysis to elucidate the genetic basis for intraserogroup O-antigen structural variations. We assembled a database of O-antigen gene islands from over 11,000 sequenced Salmonella strains, revealing (i) that gene duplication, pseudogene formation, gene deletion, and bacteriophage insertion elements occur ubiquitously across serogroups; (ii) novel serotypes in the group O:4 B2 variant, as well as an additional genotype variant for group O:4, and (iii) two novel O-antigen gene islands in understudied subspecies. We thus comprehensively defined the genetic basis for O-antigen diversity.IMPORTANCE Lipopolysaccharides are a major component of the outer membrane in Gram-negative bacteria. They are composed of a conserved lipid structure that is embedded in the outer leaflet of the outer membrane and a polysaccharide known as the O-antigen. O-antigens are highly variable in structure across strains of a species and are crucial to a bacterium's interactions with its environment. They constitute the first line of defense against both the immune system and bacteriophage infections and have been shown to mediate antimicrobial resistance. The significance of our research is in identifying the metabolic and genetic differences within and across O-antigen groups in Salmonella strains. Our effort constitutes a first step toward characterizing the O-antigen metabolic network across Gram-negative organisms and a comprehensive overview of genetic variations in Salmonella

Optimal Selection of Metabolic Fluxes for in Vivo Measurement. I. Development of Mathematical Methods

The measurement of uptake and secretion rates is often not sufficient to allow the calculation of all internal metabolic fluxes. Measurements of internal fluxes are needed and these additional measurements are used in conjunction with mass-balance equations to calculate the complete metabolic flux map. A method is presented that identifies the fluxes that should be selected for experimental measurement, and the fluxes that can be computed using the mass-balance equations. The criterion for selecting internal metabolic fluxes for measurement is that the values of the computed fluxes should have low sensitivity to experimental error in the measured fluxes. A condition number indicating the upper bound on this sensitivity, is calculated based on stoichiometry alone. The actual sensitivity is dependent on both the flux measurements and the error in flux measurements, as well as the stoichiometry. If approximate physiologic ranges of fluxes are known a realistic sensitivity can be computed. The exact sensitivity cannot be calculated since the experimental error is usually unknown. The most probable value of the actual sensitivity for a given selection of measured fluxes is estimated by selecting a large number of representative error vectors and calculating the actual sensitivity for each of these. A frequency distribution of actual sensitivities is thus obtained giving a representative range of actual sensitivities for a particular experimental situation